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We relate disentanglement and decoherence rates in a pair of three-level atoms subjected to multi-local and collective pure dephasing 
noise acting in a preferred basis. The bipartite entanglement decay rate, as bounded from above by the negativity, is found to be greater 
than or equal to the dephasing-decoherence rates characterized by the decay of off-diagonal elements in the corresponding full density 
matrix describing the system or the reduced density matrix describing either qutrit, extending previous results for qubit pairs subject to 
such noise. 

| 1 „ , 1 Introduction 

C ; 

^ , Quantum entanglement can be affected by a range of noise sources, both quantum and classical in nature. 

q-i Noise can give rise to a loss of entanglement over a broad range of quantum states. Recently, the relation- 

L ~ ~~ '■ ship between dephasing-decoherence and bipartite entanglement reduction under basis-specific classical 

! noise has been studied in two-qubit systems (HEJESIUEJEIIT] an< ^ m P a ^ rs °^ QU^nt systems (HJQlllCLl]; this 

' relationship has been studied in two-qubit systems under quantum dissipative vacuum noise as well [12] . 



£2^ . For example, it has been noted that initially entangled two-qubit systems can suffer "entanglement sud- 
\ den death," in which a two-qubit system may suddenly lose entanglement in a finite time, even though 
£NJ ■ each qubit itself maintains its quantum-coherence [3j|4]. Here, we present the first general analysis of the 
. disentanglement and dephasing of two qutrits as realized in atoms with "V"-type energy- level configura- 
' tion under a classical pure dephasing noise model, at the multi-local and collective level. We compare the 
t^T* ■ timescales of disentanglement and dephasing-decoherence, the latter timescale specifically in the basis on 
. . ! which this noise acts as in previous studies. 

In Sec. 2, we introduce our dephasing model. In Sec. 3, we examine the effects of the multi-local and 
^ . collective dephasing noise on a general state. Finally, in Sec. 4, we examine two specific classes of states 
^ | and compare their disentanglement and decoherence rates explicitly. The bipartite-entanglement decay 
■ - - ■ rate, is found to be greater than or equal to the dephasing-decoherence rates, characterized by the decay 
of off-diagonal elements in the corresponding density matrix describing the full system or the reduced 
density matrix describing either qutrit, when dephasing occurs. This result is the most general yet to be 
obtained in the comparative study of decoherence and disentanglement in two-qutrit systems. 

2 Model 

Our model consists of two-three level systems subjected to the time-dependent Hamiltonian 



H{t) = -f [bl } it) < + (t) a* + (t) (at + a?) \ , (1) 



d 1 ) U\ , , (1) M B , , (2) M / A 
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where cr z = diag(l, e~, e~) is the dephasing operator for three-level systems with subscripts denoting 
qutrits A, B, or both, h = 1, the time dependent noise terms bv{i = 1, 2) refer to statistically independent 
classical Markov processes satisfying (bx (t)) = and (bx (t) bx (t')) = j^S (t — t') , with X = A, B; 
(6ab (t)) = and (6ab (t) &ab (*')) = ~~ *0> (' ' ') * s the ensemble time average, Ti and denote the 

phase-damping rates associated with bx it) (X = A, B) and &ab (t), respectively. 

The one-qutrit standard-basis eigenstates are {|0), | + 1), [ — 1}}, for example, representing the ground 
state, first-excited state, and second-excited state of the atom, respectively. We assume that the states 
I + 1},| — 1) couple to the ground state but not to each other. Here, we notate the standard two- 
qutrit basis states via the obvious one-to-one correspondence {|1), |2), |3), |4), |5), |6), |7), |8), |9}} = 
{|00>, |0, +1), |0, -1), I + 1, 0), I + 1, +1), I + 1, -1), I - 1, 0), I - 1, +1), I - 1, -1}}, for simplicity ' 

The time-dependent density matrix for the two-qutrit system is obtained by taking ensemble av- 
erages over the three noise fields, b\(t), &b (i), ^AB (*), that is, p{t) = (Pst (*))a(b AB) j where the 
statistical density operator p s t (t) and the unitary operator U (t) associated with H (t) are p s t (t) = 
U (t) p (0) W (t) and U (t) = exp — i J* dt'H (t r ) , respectively. It is helpful to consider the dynamical 

evolution of p(t) as a completely positive trace preserving (CPTP) linear map £(p), that is, a combination 
of local and collective quantum channels, any of which can be turned off in particular cases, taking an input 
state p (0) to the output state p (t) given by the operator sum p(t) = £ (p (0)) = Y^,u=i ^\ (*) P (0) ^ (*) 1 
where are decomposition operators that satisfy the completeness relation E^E^ = I . In each 
of the various cases considered here, the internal structure of the Eu in accordance with the Hamil- 
tonian; various terms may or may not nontrivially contribute in a given case. The most general solu- 
tion of p(t), assuming that the system is not initially correlated with any of the three environments, is 

p{t) = EJj=iELi ( D k BfE f jE t f )p(°) ( E t E f D k B )> where the terms describing the interaction with 
the local magnetic fields 6a {t) and 6b (t) involve the decomposition operators E A = diag(l, 7a (t) , 7a {t))<S> 
h,E$- = diag(0,w A (t),0) ® h,Et = diag(0,0,u;A(t)) ® h,Ef = h ® diag(l )7 B (t) ,7b {t)),Ef = 
I3 (2) diag(0,a;B (t) ,0),and Ef = I 3 ® diag(0,0,ci;B (t)), and the terms associated with the global mag- 
netic field &ab (t) involve the operators D± B = diag(7 A B (t) , 1, 1, 1,7ab (t) , 1, 1, 1,7AB (t)), D AB = 
diag(w A Bi(t),0,0,0,u;AB2(t),0,0,0,cuAB2(t)), ^ B = diag(0,0,0,0,u;AB3 (t) , 0,0,0, wabs (*)) • The 
time-dependent parameters appearing in the matrices above are given by: 7, (t) = e~ t ^ 2Ti , u>i (t) = 

0^,7(*) = e-^^ait) = Jl^ti,wa(t) = and u> l3 (t) = J (l - J) (l - 7 , 4 ) , 

where Tj = ^-(i = 1,2) are the phase-relaxation times associated with the the pertinent qubits A and 
B, respectively, as introduced in [8], with Ti being the rate parameters. From here on, for tractability of 
notation, time does not explicitly appear as an argument for these quantities but is implied. 

Dephasing-decoherence rates are characterized by the decay of off-diagonal elements in the corresponding 
full density matrix describing the system or the reduced density matrix describing either qutrit given via 

characteristic decay times. Entanglement is bounded from above by the negativity J\f (p) = ^ p A ^ 1 — -, 
where p TA is the partial transpose of the density matrix with respect to qutrit A and denotes the 
trace norm |13j . 



3 General Case 

In the standard-basis representation of Eq. [21 the generic pure state of the two-qutrit system is | x I / )ab = 
ai 1 1) + a2|2) + S3 1 3) + 04 [4) + 05(5) + a6|6) + 07(7) + as\S) + ag|9), a normalized state- vector with Oj € C 
and Yli=i®i = 1- Our analysis proceeds as follows. We find the explicit time evolution of the general 
state subjected to noise from the multi-local dephasing channel £T and the collective dephasing channel 
T>. The decoherence timescales are then determined, as characterized by the decay of the off-diagonal 
elements at the level of the full density matrix of two-qutrits, as well as the reduced density matrix 
of each individual qutrit, for each of the dephasing cases, multi-local and collective. We then analyze 
disentanglement timescales, using the monotone of negativity, when there is decoherence for the general 
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state in each of these cases. In order to compare decoherence and disentanglement behavior explicitly, we 
then consider the the behavior of two particular subclasses of states, the robust class and the fragile class. 

First, consider the most general initial two-qutrit pure state, pab (0) = P(|vP)ab)j where P(|*)ab) 
is the projector corresponding to the generic state-vector argument |^)ab explicitly given above, under 
multi-local and collective dephasing noise. 



3.1 General Case: Multi-Local Dephasing Channel ST 



G,SF / .x 
Pab W 



(\a\\ aia^B aia37B aia^TA 01O57A7B oiOg7A7B oio 7 7A aiOg7A7B aiOg7A7B \ 

a 2 a\7B |a 2 1 2 02^37! 02al7A7B a 2 ailA a 2 alj A 7 B «20 7 7A7b a 2 a^7A a2a 9 7A7B 
a3ai7B a 3 a* 2 ^ I«3| 2 a3al7A7B a3a57A7i «3G67A 03a 7 7A7B a 3 ag7A7B «3«97A 
04fli7A 04027A7B 04037A7B | «4 1 2 a40s7B 04067b a 4 a 7 7 A a 4 ag7 A 7B a 4 Og7 A 7B 
a5ai7A7B a5fl27A a 5 a§7A7l a5al7B |a 5 1 2 a 5 a§7| a 5 a 7 7 A 7 B 05037! a5a 9 7A7i (j 2 ) 
06017A7B 06027A7b aea^A a6«47B 06057b l«e| 2 a 6 a77 A 7B a 6 ag7 A 7| a 6 Og7 A 
a 7 a^7A 07027A7B 07037A7B 07^171 ^Klx 7b 07067a7b (07 1 2 a 7 a|7 B 070^78 
«8al7A7B 08027A 08037A7b «8047a7b a 8 a57 A a 8 ag7 A 7| a 8 a 7 7B |a 8 | 2 a 8 a* 9 ^ 
\a 9 a^7A7B a9«27A7B «9«37A a 9 a|7 A 7 B a 9 a^7 A 7| a 9 a^7 A a 9 a 7 7 B aga^l |a 9 | 2 

is the time-evolved full density matrix for multi-local dephasing. Our density matrices are labeled from here 
on as p^f , with X denoting the pertinent qubits (A, B, AB), c is the class examined (G,F,R) as described in 
the next section, and C denoting the pertinent dephasing channel (£, T, V ) As noted above, for tractability 
of notation, time does not explicitly appear as an argument for the 7a and 75. Note that one can recover the 
local dephasing channel £{T) for qutrit A(B) by effectively freezing the time parameter in the exponential 
of the 7b (7a) factors of the opposite channel. The effect of two local dephasing noise sources, one for 
each qutrit, is simply the combination of independent effects of the individual local dephasing channel 
just described. Thus, the combined matrix is just the component-wise multiplication of exponential decay 
factors arising from each of the local dephasing channels. Due to their appearance in exponents, the decay 
rates add. 

Let us designate the different timescales as T^^ ecC and r^^, representing the decoherence and 
disentanglement timescales, respectively; i denotes the number of qutrits affected (i = 1,2), C denotes 
the pertinent dephasing channel {£,J-,T>), j is used as an index further to discriminate the timescales, 

and c denotes the class examined (G,F,R). The differing timescales of reduction of off-diagonal elements 

G( 1) 

consist of various combinations of 7a and 75. Here, we have the four timescales of decay as T 2 _^ ec£ jr = 

2(rr), ^Scc,^ = (rr). T 2-dec,£F = 2 (3r7)> and T 2-dec,£T = (srt)- The reduced density matrices of the 
individual qutrit subsystems are 

(|oi| 2 + |«2 1 2 + |a 3 1 2 (aial + a 2 a* 5 + a 3 ag)7A (aia 7 + a 2 ag + a 3 ag)7 A ,\ 
(a 4 a^ + a 5 a* 2 + a 6 a* 3 ) 7A |a 4 | 2 + |a 5 | 2 + |a 6 | 2 (a 4 a 7 + a 5 a* + o 6 o§) 7 A I 3 ) 
(a 7 a\ + a 8 a2 + 0903) 7a (a 7 a| + a 8 al + a 9 a@) 7 A |a 7 | 2 + |a 8 | 2 + |a 9 | 2 / 

(|ai| 2 + |a 4 | 2 + [a 7 | 2 (0102 + 0405 + a 7 a§) 7b (aiOg + a 4 ag + a 7 ag)7B\ 
(a 2 a* + a 5 a| + a 8 a 7 )7 B |a 2 | 2 + |a 5 | 2 + |a 8 | 2 (0^3 + 0505 + 0309)7! (4) 
q 2 2 2 / 

(a 3 a\ + a e al + a 9 a 7 ) 7 B {a 3 a 2 + a 6 al + a 9 ag) 7^ |a 3 | + |a 6 | + |a 9 | / 



Again, one finds differing dephasing decoherence timescales. In particular, the off-diagonal elements decay 
in two different timescales before the density matrix becomes fully diagonal in the basis under consideration, 
T G W - 2(-L) and r G(2) - f-M 
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The negativity, which bounds the (non-bound) entanglement, 



1 + 



ai| 2 + ([a 4 | 2 + |a 7 | 2 )7il fl«i| 2 + d«2| 2 + K^bl 



+ 



+ 



+ 



+ 



+ 



+ 



a 7 



ai 



«3 



«6 



Of) 



«2 



" S 



«5 



2 , I- i2 2 i I- |2 4 

+ |oi| 7a + l°4| 7a 



|a 7 |^ + (|a 8 |^ + |a9| 2 )72 



2 , I- |2 o t I- |2 4 

+ |ai| 7a + \ a 7\ 7a 



I- |2 , /i- |2 . I- i2\ 2 

|a-4 1 + (|o5| + |a 6 | )7r 



+ (|a 6 | + |a 9 | )t 



a 



I- [2 | [ — i2 2 i I- |2 4 

|«3| + 7b + l°2| 7b 



2 , l- |2 2 t I- |2 4 

+ |a-3 1 7a + l°9l 7a 



I- [2 | [ — i2 2 , t- |2 4 

|a 6 | + |a 4 | 7b + |a 5 | 7g 



2 , i- i2 7 i I- |2 4 

+ |o 3 | 7a + Kl 7a 



I- [2 i [ — |2 2 , I- |2 4 

|«9 1 + Kl 7b + l°8| 7b 



+ (|a 5 | + [a s ] )7 



A 



I- [2 i [ — |2 2 , t- |2 4 

|«2 1 + |ai| 7g + |a 3 | 7g 



2 , i- |2 2 i I- |2 4 

+ |o 2 | 7a + l a sl 7a 



I- |2 i I - 1 2 2 i I- |2 4 

|a-s | + \a 7 \ 7g + |a 9 | 7 B 



2 i i- |2 2 t I- |2 4 

+ |a 2 | 7a + l a sl 7a 



I- |2 i I - 1 2 2 , t- |2 4 

|o 5 | + |a 4 | 7b + \ae\ 7b 



(5) 



is rather unwieldy. It is therefore illuminating to characterize its behavior for each of the dephasing 
cases: no dephasing, one-qutrit local dephasing, and two-qutrit multi-local dephasing. In the expression 
for negativity, we can isolate the behavior of one local dephasing channel with respect to the other by 
setting the decay factors corresponding to the other channel to one, effectively freezing that factor to time 
zero. 

(1) In cases where there is no dephasing, any decay factors multiplying probability amplitudes do not fall 
off exponentially but instead their value remains at unity. J\f (p (0)) = J\f (p(t)) = 1, which corresponds to 
a maximally entangled state for all time. Thus the (trivial) dephasing and decoherence rates are equal. 

(2) In the case of a single local dephasing channel, either channel £ on qutrit A only or channel T on 
qutrit B only, the negativity J\f (p(t)) — > c > in the large-time limit. 

(3) Finally, in the case in which there is dephasing on both local channels £ and J 7 , J\f (p(t)) — > in the 
large-time limit. 

One then notes that in all these cases decoherence never proceeds faster than disentanglement, as we 
previously found to be the case for qubit pairs [7]. In later sections, when examining specific classes of 
states, we compare decoherence and disentanglement more explicitly. 
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3.2 General Case: Collective Dephasing Channel T> 

When subjected to collective dephasing noise, the time-evolved density matrix of the two-qutrit system in 
general case is given by 



Pab (*) 



/N 2 

a 3 a*7 
04 a \ 7 

— — * 4 

05^1 7 
a 6 a|7 

a 7 a^7 

a 8 a*7 

\a 9 ai7 4 



010^7 01037 

I- i2 - 
a 2 



a 3 a 2 
050^7 

O.gO.2 
0,70,2 
agOg 
090^7 



0203 

I- |2 
l a 3| 

0403 

05037 

OgOg 

a 7 a* 3 
a 8 a^ 
09O37 



aia 4 7 

a,2a| 

a 3 a| 

I- |2 
I a 4 1 

a 5 a|7 
a 6 a| 
07 a| 
a 8 a| 
a 9 a|7 



«1«57 
02057 
03057 
04O57 

I- |2 
l a 5| 

O6O57 
07057 

a 8 a^7 
0905 



aiag7 

a 2 ag 

a 3 ag 

0405 

050^7 

I- |2 
|« 6 1 

o 8 a^ 

090g7 



aia 7 7 

0207 

a 3 a 7 

0407 

a 5 a 7 7 

060 7 

I- |2 
07 



oia 8 7 

020 8 

a 3 ag 

040g 

a 5 ag7 

a 6 ag 

a 7 ag 

I- |2 
|«S I 



a 8 a 7 
a 9 a 7 7 aga 8 7 



oiOg7 4 \ 
02097 
o 3 ag7 
04097 
a 5 ag 
a 6 a^7 
07097 



a 8 a§7 

1- |2 
a 9 



(6) 



Note the form of the matrix, in which there are regions unaffected by dephasing, with decoherence effects 
along the edges. We show in the next section that this gives rise to a class of states that can exist in these de- 
coherence free subspaces, where entanglement is also be explicitly characterized. We see in the above general 

G(l) 

matrix that there exist two timescales of off-diagonal element decay, T 2 _ d g C -p 
The reduced density matrices are given by the following matrices, 



ari( 3 _ (J_ 

ana r 2 _ d v - { 2 r 2 



Pa (*) 



P B (*) 



rp G,V 



(t) 



(I- |2 1 1- |2 . 1- i2 
|oi| + \a 2 \ + |a 3 | 
a 6 ag + 040*7 + 050^7 
a 8 a2 + 070*7 + 09037 

(I- |2 1 I- |2 . I- |2 
|Ol| + |04| + |a 7 | 
a 8 a 7 + a 2 a*7 + 05047 
a 6 a4 + a 3 a*7 + a 9 a 7 7 



a 3 < + aia|7 + a 2 a*,7 

1 _ 1 2 I _ 1 2 1 - i2 
|«4 1 + |«5 1 + |«6 1 

0704 + 03057 + a 9 ag7 



0703 + 010^7 + 04057 

I- |2 1 I - |2 1 I- I? 

H- 1*3-5 1 + |a 8 | 



\ a 2\ + 1*3-5 1 + |o 8 | 0203 + 05(3 

o 3 a2 + 05057 + agag7 |a 3 | 2 + |a 6 



o 2 ag + aia 7 7 + a 3 ag7\ 
0407 + a 5 Og7 + a 6 ag7 (7) 

I- |2 1 1- |2 1 I- i2 / 
I07I + |o 8 | + |Og| / 

04 Og + 01037 + 070g7\ 

a 2 a| + 050^7 + a 8 a*,7 .(8) 

I- |2 1 I- |2 . I- |2 / 
:0 3 | + 1*3-6 1 + |<3g I / 



Note that the single-qutrit reduced density matrices do not become entirely diagonal. The reason, as we 
stated above, is the existence of decoherence free subspaces. In this case, there exists only one timescale 
in which a subset of off-diagonal elements decay, t^AccV = 

The expression for negativity for the general class under the collective dephasing channel, is similar in 
character to the negativity found in the previous subsection for the multi-local dephasing channel given 
by Eq. 5. Therefore, for clarity, we defer our discussion of the behavior of the negativity in the case of the 
collective dephasing channel until after treating specific states to the next section. 



4 Specific Classes 

In order better to distinguish decoherence and disentanglement behavior, let us now turn our attention to 
two specific classes of states. In the standard-basis representation, the generic pure state of a two-qutrit 
system is = ai|l) + 02(2) + a 3 |3) + 04)4) + 05(5} + oe|6) + 07)7) + a 8 |8) + ag|9), a normalized state 
vector with Oj £ C. The generic class of two-qutrit pure states represented by |^) contains two subclasses 
of interest, distinguished by their coherence behavior in large timescales under collective dephasing noise, 
that is, dephasing in which each qutrit interacts with the same collective noise &ab (*)• O ne class is seen 
to be fragile, whereas the other is robust. 

(i) The fragile class \(f>) = oi[l) + a5|5) + ag|9), in which 04,05, and a 9 may be non-zero and all other terms 
en = 0, has the forms 
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\<j> l ) = a^l) + a 5 |5) , (9) 
|02) = oi|l) + afl|9) , (10) 
\<fo) = a 5 |5)+a 9 |9) . (11) 

(ii) The robust class \ip) = a,2\2) + 03(3} + 04(4} + o,q\Q) + aj\7) +as|8), in which all Hi listed may be non-zero 
and a\ = ag = 0, has the forms 

= a 2 |2) + a 4 |4) , (12) 
fofe) = a 3 |3) + a 7 |7) , (13) 
|^ 3 > = a 6 |6) + a 8 |8) . (14) 

The specific forms of the two classes above constitute Bell-like states in the Hilbert space of our two-qutrit 
system, similarly to the qubit-pair-state classification given in [8]. 



4.1 Fragile Class 

The initial density matrix representing the two-qutrit fragile class is given by 



/ \ai\ 2 ■ ■ ■ aia^ • ■ ■ aiag \ 



p F AB (O) = P(|0» 



05^1 



\a g al 



l fl 5| 



a 9 a 5 



a 5 a g 



I- |2 / 



(15) 



where the dots denote the remainder of the density matrix, which is filled out by zero entries. Let us 
turn our attention to the behavior of this state when subjected to the multi-local dephasing 8 J- acting on 
qutrits A and B. Recall that the multi- local dephasing channel captures the features of each of the local 
dephasing channels fully, so for this analysis, we consider the multi-local dephasing channel only, which 
provides information about each of the individual dephasing channels if one simply turns off one of the 
channels. 



4.1.1 Fragile Class: Multi-Local Dephasing Channel £T ' . The density matrix subjected to noise is 



Pab (*) 



l a i 



a5al7A7B 



\a 9 a^7A7B 



aia^7A7B ■ • • a 1 a^ A 7B\ 



\ a 5\ 



°9 a 57A7| 



05^97171 



I — I z / 



(16) 



which has two different decoherence times of T^_ dec£: p = (jr-) and T~2_ dcc£:F = (217)- By compari- 
son, the disentanglement timescales deriving from the negativity, jV (p^B (*)) = ^ ^ 1 = ^ P -5 1 — - = 
(|ai[ |a 5 | + |ai| \a 9 \) 7a7b + |as| l«9 1 7a7 b are r Jii^ = (17) and T dii 2 ]^ = (217)- 
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Produced (t) = Tr B PAB^ (*) = Tr APAB^ (*) 




(17) 



We see that they are both initially fully mixed, so that no further decoherence is possible. 

Comparing timescales we see that decoherence never proceeds faster than disentanglement, r^ is £:F < 

T 2-dec,ET and T dis£F - T l-dec,£T 

4.1.2 Fragile Class: Collective Dephasing Channel T>. The two-qutrit density matrix, as given by 



Pab (*) 



l a i| 



— — * 4 



\a 9 a^7 4 



aia* 5 j 4 ■ ■ ■ aia^7 4 \ 



l a 5| 



a 9 al 



aaa 9 



I- |2 / 

\a 9 \ j 



(18) 



decays according to a single timescale, T^decX" = (^r)- The disentanglement, characterized via the nega- 
tivity, M (p^B (*)) = ~ it — " = ~ it — " = 1^1 I'M + (l^il 1^1 + l^il l^9l)7 4 - Disentanglement proceeds 
at only a single timescale, rJ isV = (^)- 
The single qutrit matrices are always 



Produced (t) = TtbPab (*) = Tv APab (*) 



ai| 







Ids I 







(19) 



|a 9 | 



Again, we see that the single qutrit states are both fully mixed, so that no further decoherence is 

possible. Comparing timescales, we see that decoherence never proceeds faster than disentanglement: 
F ^ F 

T dis,X> — T 2-dec,X>- 



4.2 Robust Class 

The two-qutrit density matrix for the robust class, under the multi-local dephasing channel, is given by 

/0 000 o o\ 

|a2| 2 02^3 0204 a2(lg a2a 7 a2a§ 
a^a^ 1 03 1 2 0304 a^al 0,30,7 0,30% 

(24(22 (14(23 1 0,4 1 2 (240,g (240.7 (24 (2g 

Pab (0) = . (20) 

______ I _ 1 2 — — — — 

a§a,2 oqo^ a§a\ |ae| «6°7 «6 a s 

(27(22 07(23 0704 a>ag |(2>| a7(2g 

2 

asaJj a8«3 fl8«4 agag a%ctj \a%\ 



yo o 







0/ 



The rows and columns of zeros in this matrix are exactly the rows and columns that are affected by the 
collective dephasing channel, giving rise to the robust character of this density matrix. For analysis in 
this section, we consider, without loss of generality, the first form of this class = 02 12) + 04 14), which 
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Jaeger and Ann 



allows the decoherence and disentanglement timescales to be seen clearly. If we were to include the entire 
subclass, we would have the same behavior as our forthcoming analysis. 



4.2.1 Robust Class: Multi-Local Dephasing Channel £T ' . 

multi-local dephasing channel acting on is given by 



The time-evolved density matrix for the 



R,£F 

Pab 



(t) 



/0 








0\ 


|a 2 | 2 


a 2 a^7A7B 











00 


a 4 a^7A7B 


1- |2 
a 4 





















































\0 








0/ 



(21) 



We see that the off-diagonal elements decay according to the timescale tJ^_ 
glement, characterized via the negativity, is given by, TV (pW (t)) 



2-dec,£J r 

-l Wp 1 



(jt|-). The disentan- 



«2 O4 7A7B- 



Disentanglement proceeds at the same single timescale, r^ s £:F = (j^)- The re duced density matrices are 
given by 



R,£F 
Pa 



(t) 



Tr B p AB 



I- |2 n 
02 



la.4 1 2 



(22) 



Pb v) 



IWaTT (t) I 




I- I 2 n 
O4 (J 



|a 2 | 2 



(23) 



always fully decohered as in the case of the fragile states. Again, we see that decoherence never proceeds 
faster than disentanglement, so that clearly the inequality r^ s v < T^_ dec v is satisfied. 



4.2.2 Robust Class: Collective Dephasing Channel T>. Here, p(t) = p (0) as defined above in Eq. 
[20l and the state neither decoheres nor disentangles: the off-diagonal elements remain unreduced and the 
negativity remains at the maximum value of 1 for all times. These states are robust against collective 
dephasing noise. 



4.3 Conclusions 

We have shown for this model of two-qutrit systems under the chosen basis-specific dephasing noise, 
whether local, multi-local or collective in nature, that disentanglement proceeds at least as fast as deco- 
herence. This result accords with previous studies of similar nature for qubit systems. Because the case of 
two-dimensions in quantum mechanics is often a special one in quantum mechanics, for example, in the 
cases of the foundational theorems of Gleason and Kochen and Specker, it is always important to examine 
cases beyond that of qubits for the possibility of different behavior from that exhibited in the special 
case of qubits. Our results suggest that the relation between basis-dependent dephasing decoherence and 
disentanglement found here can be expected to hold for pairs of initially entangled qu-d-it systems under 
analogous dephasing noise for general values of d. It remains an open question whether or not this phe- 
nomenon can be generalized to an arbitrary multi-partite quantum systems, each of arbitrary dimension, 
however, not least of all because good entanglement measures for mixed states of such systems is lacking. 



Disentanglement and decoherence in a pair of qutrits under dephasing noise 
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